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ABSTRACT : The durability and damage tolerance of laminated composites are critical 
design considerations for airframe composite structures. Therefore, the ability to model 
damage initiation and growth and predict the life of laminated composites is necessary to 
achieve structurally efficient and economical designs. The purpose of this research is to 
experimentally verify the application of a continuum damage model to predict progressive 
damage development in a toughened material system. Damage due to monotonic and 
tension-tension fatigue was documented for 1M7/5260 graphite/bismaleimide laminates. 
Crack density and delamination surface area were used to calculate matrix cracking and 
delamination internal state variables to predict stiffness loss in unnotched laminates. A 
damage dependent finite element code predicted the stiffness loss for notched laminates 
with good agreement to experimental data. It was concluded that the continuum damage 
model can adequately predict matrix damage progression in notched and unnotched lami- 
nates as a function of loading history and laminate stacking sequence. 

KEY WORDS: composites, graphite/epoxy, damage, matrix cracks, delamination, inter- 
nal state variables. 
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INTRODUCTION 

B ecause nr thhir light weight and high specific stiffness, laminated con- 
tinuous fiber-reinforced composite materials are being used in some primary 
load bearing components in aircraft structures. However, when subjected to high 
service loads, environmental attack, impact, or a combination of any or ail of the 
above, laminated composite materials develop damage. As the service load or the 
time in service increases, damage develops and grows, becoming more severe, 
and could eventually lead to catastrophic failure. 

There are four main types of damage. These are matrix cracking, fiber-matrix 
interface debonding, delamination, and fiber fracture. Usually, matrix cracking 
and fiber-matrix interface debonding are the first forms of damage to occur, fol- 
lowed by delamination, and finally fiber fracture resulting in catastrophic failure. 
While matrix cracking is usually arrested at the fibers or adjacent plies, it will 
result in a redistribution of load to the surrounding regions. As a result, these 
surrounding regions contain stress fields which are favorable to the initiation and 
propagation of additional damage. During the accumulation of subcritical dam- 
age, changes in material stiffness and strength results in the load redistribution 
until the principle load-carrying plies are unable to support the load, in which 
case, catastrophic failure occurs. 

The initiation and propagation of damage is one of the problems in using lami- 
nated continuous fiber composite structures. To address durability and damage 
tolerance requirements, damage must be modelled and methods developed to pre- 
dict the residual strength and life of composite structures. For example, one of the 
most complicated structural configurations is that of built-up laminated compos- 
ite structures connected by mechanical fasteners such as rivets. These laminates 
with fastener holes develop local damage that cannot be easily treated using stress 
concentration factors. Another example is the non-visible damage that develops 
during foreign object impacts and ground handling accidents. Current methods 
for treating these local damage details are empirical and very conservative. 
Therefore, an accurate model of the damage initiation and propagation will en- 
hance current analysis and design capabilities thus leading to improvements in 
structural efficiency. 

Many damage progression models are being developed to model damage and 
predict life. An example of the type of microcrack damage that is currently being 
studied by damage models is shown by the X-ray radiograph in Figure I. This 
damage is both stacking sequence dependent and loading history dependent. An 
overview of damage resulting from fatigue loading in composites has been pre- 
sented by Reifsnider [1,2]. Some researchers have tried to model this damage by 
considering each crack as an internal boundary and the stress or displacement 
fields are obtained either in closed form or numerically, such as in finite ele- 
ments. This approach works well as long as there are a relatively small number 
of cracks. Talug and Reifsnider [3] have obtained finite difference approximate 
solutions to equilibrium equations to solve for interlaminar stresses in composite 
laminates. The "damage tolerance/fail safety methodology” developed by O’Brien 
[4] is an engineering approach to ensuring adequate durability and damage toler- 
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ancc by treating only delamination onset and the subsequent delamination accu- 
mulation through the laminate thickness. Chang [5 1 developed a progressive 
damage model for notched composite laminates subjected to monotonic tensile 
loading. This particular model assesses the damage and predicts the ultimate ten 
sile strength in laminates with arbitrary ply-orientations via an iterative combina 
lion of stress analysis and failure analysis. Chamis |6] studied structural charac- 
teristics such as natural frequencies and buckling loads and the corresponding 
mode shapes during progressive fracture of angle-plied polymer matrix compos- 
ites. This study concluded that the individual nature of the structural change was 
dependent on laminate configuration, fiber orientation, and the boundary condi- 
tions. The model proposed by Talreja [7-9] incorporates internal state variables 
(ISVs) for matrix cracks and delaminations and exhibits ply stacking sequence 
dependence. The ISVs are strain-like quantities which represent the damage as 
volume averaged quantities, i.e., a continuous medium. 

The treatment of a damaged volume of material as a continuous medium and 
the representation of the damage with averaged quantities was first proposed by 
Kachanov [10] in 1958 and is referred to as continuum damage mechanics. From 
continuum damage mechanics has evolved the continuum damage model devel- 
oped by Allen and Harris 1 11 — 16). This damage model utilizes ISVs and is phe- 
nomenological; however, it is formulated at the ply and sublaminate level and 
therefore accounts for the influence of stacking sequence and geometry. 

The goal of the research presented in this paper is to apply the Allen and Harris 
model to a toughened matrix composite system and to experimentally verify the 
current predictive capability of the model. The IM7/5260 material system was 
selected for characterization because the 5260 bismaleimide is a candidate 
elevated temperature polymer for the next generation supersonic transport air- 
plane. Both cross-ply and quasi-isotropic laminates were investigated in un- 
notched and open hole specimen configurations. Tension-tension fatigue loading 
was used to initiate and propagate damage at several different constant amplitude 
stress levels. No new mathematical models were developed for this study. This 
investigation was intended to be an interrogation of the suitability of the 
mechanics framework and the current capabilities of the model to predict damage 
growth in laminated composites. 

THE ALLEN AND HARRIS MODEL 

The damage model of Allen and Harris [11-16J was developed to model the be- 
havior of microcrack damage in a brittle epoxy material system by predicting 
stiffness loss and damage-dependent ply level stresses in a laminate. A summary 
of the model can be found in the literature [17 J. The model, which neglects edge 
effects, uses internal state variables to represent the local deformation effects of 
various modes of damage. Loading history dependence is modelled by ISV dam- 
age growth laws. The progression of damage is predicted by an iterative and in- 
cremental procedure outlined in the flowchart shown in Figure 2. This entire pro- 
gressive failure analysis scheme has been implemented into the NASA 
Computational Mechanics Testbed (COMET) [18]. The first block of Figure 2 is 



284 


Timothy W. Coats and Charlrs E. Harris 



Figure 2. Progressive failure analysis scheme. 
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a description of the information needed as model input. A FORTRAN code con- 
sisting of the damage dependent constitutive model and a damage growth law for 
matrix cracking was incorporated into a classical lamination theory analysis to 
produce effective lamina and laminate properties for unnotched laminates. The 
program is called FLAMSTR (Fatigue LAMinate STRess) 1 19 J and makes up the 
first constitutive module. The fourth block is a damage dependent finite element 
analysis code (18] from which the second constitutive module performs a ply 
level elemental stress analysis and simulates damage growth via damage growth 
laws for each element. Updating the damage state, block seven, for notched lami- 
nates is the result of the damage evolution calculations in block six. For un- 
notched laminates, only the first constitutive model is needed to update the 
damage state. 

MATHEMATICAL FRAMEWORK FOR DAMAGE 
DEPENDENT LAMINATION MODEL 

The damage model represents each mode of microcrack damage by the volume 
averaged dyadic product of the crack face displacement, n, , and the crack face 
normal, n,, as originally defined by Vakulenko and Kachanov (10), 



where a, y is the second order tensor internal state variable, ,v c is the crack surface 
area, and V L is the local representative volume, i.e., all stresses, strains, and ISVs 
are averaged over a local volume element. This product can be interpreted as ad- 
ditional strains incurred by the material as a result of the internal damage. The 
local ply level damage-dependent stress-strain relationship is given by |I2|, 

= Q-jH l f »! _ <*»,)/ (2) 

where o iJl are the locally averaged components of stress, Q ljk , is the ply level 
transformed stiffness matrix, t 4/( are the locally averaged components of strain, 
and a kll are the components of the internal state variable. From here on out, the 
locally averaged subscript, L, will be dropped and all quantities are considered 
locally averaged unless specified otherwise. This ply level stress-strain relation- 
ship can be implemented directly into lamination theory. The laminate equations 
are constructed by assuming that the Kirchoff-Love hypothesis may be modified 
to include the effects of jump displacements i", and u" as well as jump rota- 
tions (3 j and tp® for the / th delaminated interface. Therefore, the displacement 

field in the x, y, and z directions is defined as |I2| 

ulx.y.z) = u°(x,y) - z\0° + H(z - z,)ii l ’\ + H(z - z,)u 

rix.y.z) = v°(.v, y) - z[p° + H(z - z<W‘ I + H(z - z,)y" 

u ( v, y, c) = w"(x,y) + //(~ - j, )m " 


(3) 
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where the superscripts o imply undamaged midsurface quantities, and H(z — z.) 
is the Heavyside step function. Locally averaging the displacement field provides 
the means to obtain the locally averaged components of strain as |12,20| 
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which can be substituted back into the ply level damage-dependent stress-strain 
relationship [Equation (2)]. The laminate midplane forces and moments per unit 
width of region V L in the laminate are given by [20] 
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Therefore, the resulting laminate equations are [12] 


\n\ = £ \QKz k - z*.,)k i - \ E \QU(zi - z* 2 -,)ui + 
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- - cf-.)l« M U (7) 

4 - I 

where \Q\y is the undamaged transformed modulus of the Ath ply in the lami- 
nate, is the ^-coordinate of the Ath ply interface, e° are the locally averaged 
midsurface strains, and x are the locally averaged midsurface rotations. Finally, 
are the averaged transformed stiffness matrices related to the / th delami- 
nated ply interface, and the superscripts D and M represent delamination and 
matrix cracking, respectively. 

From Equation (I), there are nine components of the internal state variable for 
any mode of damage. However, most of these components arc zero. For intraply 
matrix cracks the resulting non-zero damage variables are <* 22 , oft, and a". 
Since each ply has in plane symmetry, a 32 is negligible compared to the other 
components [ 12,21 J Delamination internal state variables are also defined by 
volume averaged displacements. Interface displacements, w'\ r'\ and w ? at the 
i th delamination sites are represented by a?,, and aft, respectively. Delami- 
nation induced rotations about the r and y axes, yf? and 0”, are represented by 
ol 4 , and of", respectively. All other components are zero [12,21], 

Mode I Matrix Crack Growth Law 

A damage growth law has been developed for mode I matrix cracks, oj', 
where the displacement of the crack face is in a direction parallel to the crack 
face normal, i.e., perpendicular to the plane formed by the ply. The mode I strain 
energy release rate is calculated from the following equation 1 19,21], 



( 8 ) 
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where dn&fds is a linear function of the local ply stresses and describes the 
change in the internal state variable for a certain change in crack surface area, 
and can be written as 


daZ 

8x 


(</para)(a 22 ) 


(9) 


where </para is a material dependent parameter 1 19) which must be determined 
experimentally. Therefore, the damage growth law for uniaxial cyclic loading is 
given by 1 16,19) 


dN 


daV 2 

ds 


kG\ 


( 10 ) 


where k and r/ are experimentally determined material dependent parameters, N 
is the number of load cycles, and G, is the mode I strain energy release rate for 
the ply of interest. 

Mode I Delamination Empirical Formulation 

The following relationship [12] provides an empirical formula with which to 
calculate the internal state variable due to delamination 

dal _ „ * (E™ ~ E*) (S d \ 

de, " 2Q\, + Q^\S) 

where E xo is the undamaged experimental stiffness and E* is the moduli of the 
sublaminates formed by delamination [22,23]. Furthermore, n is the total 
number of plies, S D is the delamination surface area, S is the total surface area, 
and Q f, and Qu are the transformed stiffnesses of the plies above and below the 
dclamination formed by regions of delamination at the midplane. 


MODELS FOR DAMAGE GROWTH 

The material IM7/5260 was chosen to experimentally verify the continuum 
damage model. Unlike the brittle material, AS4/3502 used to develop the model, 
IM7/5260 has a toughened resin, and it was thought that this material would test 
the applicability of the model to a toughened matrix composite. The experimental 
procedure was designed based upon the data needed to successfully use the 
model. 

Three material parameters, d para, k, and 77 , must be determined experimen- 
tally. The calculations involved in determining the model parameters have been 
previously published in detail 1 19,24] . The material parameters can be shown to 
be determined from the mode I internal state variable, a 22 . This was accom- 
plished by subjecting |0/90 2 /0|, laminates to tension-tension fatigue [24] with 
maximum stress levels at 30% ultimate stress. Edge replicas and X-rays were 
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used to measure crack surface area and determine crack density in the 90 degree 
plies. The ply level transverse and shear moduli were determined from unidirec- 
tional and |45/-45| 2 , laminates [25]. The internal state variable, a 22 . was 
calculated for each crack density measurement utilizing the following equations 
[15,19| 


a 22 


Q_ 

27 


64 £ 


( 12 ) 


t =E E l - JjT, 03) 

m " C 12t A2m - 1) 2 (2 n - 1) J + C l2lJ (y) (2/i - 1 )* 

where q is the force per unit length applied transverse to the 90 degree fibers, 2F 
and 1/25 are the layer thickness and crack density, respectively, and qI2T is 
basically the local stress (a 2 ) in the 90 degree plies. C 2222 is the transverse modu- 
lus, C, 212 is the shear modulus, and m and n are the number of iterations. Fur- 
thermore, da 22 lds is obtained by plotting a 22 as a function of crack surface area 
(Figure 3). This results in the following equations [24J 

a 22 = -1.57 x 10 6 + (2.91 x 10 ~ 4 )(s) + (1.22 X I0- S )(s 2 ) (14) 

^ = 2.91 x 10-" + (2.45 x l0- 5 )(i) (15) 


where 5 is the crack surface area. Referring to Equation (9), the first material pa- 
rameter, d para, was found as the slope of da 2 Jds vs. the local stress in the 90 
degree plies, o 2 . 

Plotting a 22 as a function of the number of cycles (Figure 4) resulted in the fol- 
lowing equations [24] 

a 22 = -4.81 x I0 6 4- (3.49 x 10“ 9 )/V - (4.77 x 10 ~*)N 2 (16) 

= 3.49 x 10‘ 9 - (9.54 x 10 ~")N (17) 


Furthermore, the mode I ply level strain energy release rate, G, could be 
calculated as 


o, - (r)(^)c.) 


(18) 
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Crack Surface Area (s), In 2 


Figure 3. Internal state variable for matrix cracks in IM7/5260 as a function of crack surface 
area. 
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Figure 4. Internal state variable for matrix cracks in IM7/5260 as a function of fatigue 
cycles. 
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since da 12 lds and o 2 are known over the given number of cycles. The ply thick- 
ness, /, is actually the thickness of the consecutive ninety degree plies, i.e., / is 
equal to two times the thickness of one ply for a [0/90 2 /0|, laminate, whereas t 
is equal to the thickness of one ply in a (0/90/0], laminate. Finally, rearranging 
the damage growth law (Equation (I0)| to the form 


ds dotiJdN 

— — L(Zy 

dN da 21 lds ' 


(19) 


from Equations (15) and (17) (da 2 ildN)l{da 12 lds) was plotted as a function of the 
mode I strain energy release rate, G,. A power curve was fit to the plotted data 
and, as a result, the final material parameters, k and 17 were obtained. The param- 
eters were then used in a fatigue laminate stress program (FLAMSTR) [19] utiliz- 
ing the damage-dependent lamination theory to predict reductions in stiffness due 
only to mode 1 matrix cracking for unnotched composite laminates. 

Other contributors to stiffness loss are delamination damage and matrix crack- 
ing due to shear. The following general equation can be used to predict reduc- 
tions in stiffness for any laminate ( 21 ] 


E M 


E 


to 


= 1 


AE M AE n AE S 



( 20 ) 


where AE X \ A E n % and AE S refer to the change in stiffness due to mode I matrix 
cracking, delamination, and matrix cracking in the 45 degree plies, respectively. 
However, we assume the stiffness loss due to the 45 degree ply matrix cracks is 
negligible and therefore ignore it primarily because previous work |24] showed 
the effects of these 45 degree ply cracks to be small compared to other modes of 
damage. Furthermore, based on results from a separate study (26), plastic strain 
in these laminates was negligible at the applied stress levels. Delamination, on 
the other hand has a dominant effect. We define A E n /E xo by the following equa- 
tion for any number of delamination sites as [ 21 ] 


AE? 

E xn 


1 

TE U 


E I Q>,M 


daZ 

Be, 


( 21 ) 


where 


, Qn + Q f, 


( 22 ) 


which, by substituting Equation (11) into Equation (21), can be reduced to 



E x 


( 23 ) 
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The use of Equation (23), which is identical to O’Briens equation for stiffness of 
a partially delaminated laminate [22], requires the experimental determination of 
the delamination area S D . As such, the model is only correlative for delamination 
damage. Furthermore, the computation of E* must be made on the basis of 
whether the delamination is an edge delamination versus a local delamination 
that initiates at a matrix crack [27]. 

MATERIALS AND EXPERIMENTAL PROCEDURE 

As previously mentioned, the material chosen to experimentally verify the con- 
tinuum damage model was IM7/5260 graphite/bismaleimide laminates. Cross- 
ply and quasi-isotropic laminates were cut into 2.54 cm X 25.4 cm (1 " X 10") 
coupons. The layups used were [0/90 2 /0],, [0/90-,],, [0/45/ — 45/90],, and 
[90/ -45/45/0], notched and unnotched. The notched laminates had a 6.35 mm 
(1/4 inch) hole drilled in the center. Each laminate was subjected to tension- 
tension fatigue up to 100,000 cycles at a frequency of 5 Hz and a stress ratio of 
0.1. Six specimens of each layup were fatigued such that three specimens were fa- 
tigued at maximum load levels of 30% of ultimate failure load and three more at 
65%. In situ edge replicas and X-ray radiographs to characterize damage were 
taken throughout the testing so that the specimen did not have to be removed from 
the load frame. A Direct Current Displacement Transducer (DCDT) with a 4" 
gage length was secured to the specimen and remained secured and untouched for 
the entire test. The fatigue loading was stopped periodically to collect edge rep- 
licas and X-rays, and to monotonically load the specimen for strain measure- 
ments. 


EXPERIMENTAL CHARACTERIZATION 

To experimentally verify the model, experimental stiffness loss directly cor- 
relating to damage accumulation was compared to model predictions. The fol- 
lowing procedure was used to document stiffness loss and characterize damage. 
The function of the DCDT was to measure increases in strain due to damage ac- 
cumulation over a 4" gage length. Initial strain measurements as well as edge 
replicas and X-rays were taken to ensure no initial damage and establish the ini- 
tial undamaged stiffness. Periodic strain measurements, edge replicas, and X-rays 
showed a degrading trend in stiffness with damage accumulation. The in situ 
edge replicas and X-ray radiographs provided a direct and accurate correlation 
between damage accumulation and stiffness loss. The edge replicas and X-ray 
radiographs provided the means to measure matrix crack and delamination sur- 
face areas. At the end of the test, the specimens were monotonically loaded to 
failure. Comparing these failure loads with the failure loads of the undamaged 
specimens provides the residual strength at 100,000 cycles for a given fatigue 
stress level. 


EXPERIMENTAL VERIFICATION OF MODEL 

The material parameters, dpara, k, and r/, are determined from the mode I 
matrix crack surface area. The unnotched [0/90 2 /0], laminates provided this in- 
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formation since the matrix cracks in the 90 degree plies are strictly mode I and 
no cracks accumulate in the 0 degree plies. The crack surface area as a function 
of fatigue cycles was measured from the edge replicas and the X-ray radiographs 
[24]. From there, the matrix crack internal state variables were calculated and 
used to determine the material parameters as mentioned previously in the section 
“Models for Damage Growth The parameters are then used in the mode I matrix 
crack growth law to compare stiffness loss predictions due only to mode I matrix 
cracking to the experimental stiffness loss of the other notched and unnotched 
laminates. For those laminates with delaminations, the stiffness degradation was 
greater. The deiamination surface area and delamination locations, as determined 
from the X-ray radiographs and the edge replicas, respectively, were used in con- 
junction with Equation (24) to provide a prediction of stiffness loss due to 
delamination. 


RESULTS 

The analysis of the experimental data provided the following material parame- 
ters: 


dp ara = 8.8686 x I0~ 7 /Ib 



7 ) = 5.5109 


These parameters were then used, via the mode I matrix cracking growth law, in 
a fatigue laminate stress program (FLAMSTR) [19] and a damage dependent 
finite element code installed in the COMET [18] to predict reductions in stiffness 
due only to mode I matrix cracking illustrated in Figure 5 for [0/90^], laminates. 
These predictions are illustrated in Figure 6 where we see that the predicted re- 
ductions in stiffness are in close agreement with the experimental stiffness loss 
due solely to mode I matrix cracking in the 90 degree plies of the [0/90,], lami- 
nates. Not only are these two distinct trends predicted from two different constant 
amplitude maximum stress levels, but the predicted trends follow the experimen- 
tal data very well. The X-ray radiographs in Figure 5 exactly correlate with one 
data point in Figure 6 for each constant amplitude maximum stress level, and two 
separate trends in stiffness loss were observed for two separate damage states. 

If delamination is also considered, damage-dependent empirical formulas may 
be used to determine the predicted reductions in stiffness due to delamination. 
The [0/45/ -45/90), laminates suffered edge delamination primarily at the 
-45/90 interface. A comparison of the delamination surface area with measured 
stiffness loss at 100,000 cycles is illustrated in Figure 7. The experimentally 
measured delamination surface areas were used in Equation (24) to predict 
stiffness loss due to delamination growth. The experimental results with their 
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Figure 5. Loading history effect on mode I matrix cracking in [0/90 3 U IM7/5260 at 700,000 
cycles. 
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predictions for laminates subjected to 50% and 60% of ultimate load are shown 
in Figures 8 and 9. A graphical comparison of these two loading histories is 
shown in Figure 10. Once again, the model predicted greater stiffness loss for a 
greater maximum stress level. Furthermore, as was expected, stiffness loss in 
these laminates was not greatly affected by the mode I matrix cracking alone. 

To further verify the delamination empirical formula for stiffness loss. Equa- 
tion (24) was applied to a different stacking sequence. A comparison of the dif- 
ferent stacking sequence is illustrated in Figure 11 where the arrows indicate loca- 
tions of primary local and edge delaminations for the [90/ —45/45/0), and 
[0/45/ — 45/90), layups, respectively. The X-ray radiographs in Figure 11 show a 
distinct difference in the damage state and a larger reduction in stiffness per 
delamination surface area for the (90/ —45/45/0), laminate. Figure 12 is the ex- 
perimental and predicted stiffness loss of a [90/ — 45/45/0], laminate. Close atten- 
tion should be paid here to the difference between £* for edge delamination vs. 
local delamination [27|. 

The analysis of the notched laminates yielded good results as well. An illustra- 
tion of the damage commonly seen in a notched laminate for two different stack- 
ing sequences is shown in Figure 13. Initially, one would expect the matrix crack 
in the zero degree ply adjacent to the hole (axial split) to have a significant effect 
in reducing the stress concentration at the notch. This would increase the global 
displacements of the laminate prior to failure. A damage dependent finite element 
code [18| within the NASA COMET was used to calculate mode I matrix crack- 
ing damage variables, laminate stresses and strains, and far field displacements. 
The finite element mesh in Figure 14 is a quarter section of the notched laminate. 
The finite element code predicts the damage state in each element as a function 
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Figure 8. Experimental vs. analytical stiffness degradation of [0/45/ - 45/90 ) s IM7/5260 for 
tension-tension fatigue at P max = 50% Pult. 
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Figure 9. Experimental vs. analytical stiffness degradation of [0/45/- 45/90], IM7/5260 for 
tension-tension fatigue at P ma * = 60% Pult. 
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Figure 10. Experimental vs. analytical stiffness degradation of [0/45/ - 45/90], IM7/5260 
for tension-tension fatigue. 
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Figure 1 1. X-ray radiograph of quasi-isotropic IM7/5260 laminates. 
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Figure 12. Stiffness degradation of [90/ -45/45/Oh IM7/5260. 
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Figure 13. X-ray radiograph of I M 7/5260 notched laminates. 
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Figure 14. Finite element mesh for laminates with a central circular notch 
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Figure 15. Stiffness loss of IM7/5260 laminates with a central circular notch. 
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of the local element stresses. An iterative procedure is used to calculate the ele- 
ment damage-dependent properties and associated load redistribution throughout 
the finite element model. The analytical far field displacements calculated over a 
4" gage length are compared to the experimental stiffness loss in Figure 15. This 
figure illustrates the ability of the code to predict separate trends in stiffness re- 
ductions due to mode I matrix cracking for different constant amplitude stress 
levels and layups for a spatially varying damage state. Edge and local delamina- 
tions were also figured into the predictions and were found to have very little ef- 
fect, especially compared to the mode I matrix cracking of the axial split for the 
|0/9(U laminate. This is attributed to the small delamination surface areas as 
well as the fact that the difference between E * and E xo is small. The comparison 
of the two laminates given in Figures 13 and 15 confirms the ability of this model 
to predict damage growth as a function of the laminate stacking sequence. The 
1 0/90,), laminate has more severe axial splitting, i.e., more mode I matrix crack- 
ing, thus the predicted loss in stiffness due to mode I matrix cracking is larger. 
The reductions in stiffness are greater for this laminate because more load is 
transferred away from the stress concentration at the hole. 

CONCLUSIONS 

The damage model originally developed by Allen and Harris was applied to the 
IM7/5260 toughened matrix material system. Experimental verification of the 
model was established by comparing the stiffness loss of cross-ply and quasi- 
isotropic laminates with and without open holes for tension fatigue loadings. The 
model has predictive capability for intraply matrix cracks and correlative capa- 
bility for delaminations. The model successfully predicted both the effects of 
laminate stacking sequence and loading history on damage growth and stiffness 
loss. The ability of the model to predict damage growth in the open hole speci- 
mens was particularly encouraging. This illustrates that the model is appropriate 
for spatially varying damage and not confined to uniform damage that typically 
develops in the gage length of an unnotched uniaxial test specimen. The spatial 
variation in damage is treated through the finite element discretization since the 
damage is assumed to be uniform within an element. The empirical formulae for 
delamination provided trends in stiffness loss that agreed with the experimental 
trends. It should be noted here, however, the predictive capability of this model 
would increase dramatically if delamination growth laws were available rather 
than the empirical relations which need the extent of the delamination to be 
specified rather than predicted. 

NOMENCLATURE 

a two times crack density 
C tJkl ply level stiffness tensor 
Jpara material parameter — slope of dotfds vs. local ply 
stress 

E* axial stiffness of a laminate completely delaminated 
along one or more interfaces 
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E xn undamaged longitudinal stillness 
EjE xo longitudinal stiffness loss 

G i mode I strain energy release rate 
(t/ l2 ),r,, experimental shear modulus at final cycle 
G \ 2 i, undamaged shear modulus 
k growth law material parameter 
n number of plies 
n\ crack face normal 
\Q] ply level transformed stiffness matrix 

average transformed stiffness matrices of sub- 
laminates formed by delaminations 
()}, stiffness at the top sublaminate at a delamination in- 
terface 

Qu stiffness of the bottom sublaminate at a delamination 
interface 

.v r crack surface area 
S (> delamination surface area 
S total surface area 
7 one-half ply thickness 
/ thickness of consecutive 90 degree plies 
T laminate thickness 
f, twice the ply thickness 
crack face displacement 
n, r, n locally averaged displacements 
u ° (.v . v) , r ° (. v , v ) , (.r , y) undamaged midsurface displacements 

n*" delamination interface displacements 
V r local representative volume 
Zk Z coordinate of the kih ply interface 
a™ matrix cracking internal stale variable 
delamination internal state variable 
delamination induced rotations 
AE x1 /E to mode I matrix cracking stiffness reductions 
A E s !E xo 45 degree ply matrix cracking stiffness reductions 
AE n IE x „ delamination stiffness reductions 

( locally averaged components of strain 
x locally averaged midsurface rotations 
r} growth law material parameter 
q force per unit length applied transverse to the 90 
degree plies 

a locally averaged components of stress 
o 2 local stress in the 90 degree plies 
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